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Thermal noise can destroy topological insulators (TI). However we demonstrate how TIs can be 
made stable in dissipative systems. To that aim, we introduce the notion of a band Liouvillian as 
the dissipative counterpart of a band Hamiltonian, and show a method to evaluate the topological 
order of its steady state. This is based on a generalization of the Chern number valid for general 
mixed states (referred as density matrix Chern value), which witnesses topological order in a system 
coupled to external noise. Additionally, we study its relation with the electrical conductivity at 
finite temperature, which is not topologically invariant. Nonetheless, the density matrix Chern 
value represents the part of the conductivity which is topological due to the presence of quantum 
mixed edge states at finite temperature. To make our formalism concrete, we apply these concepts 
to the two-dimensional Haldane model in the presence of thermal dissipation, but our results hold 
for arbitrary dimensions and density matrices. 
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Introduction. — In a recent work [1^, we have shown that 
certain one-dimensional topological insulators (TI) lose 
the topological protection of their edge states when they 
are coupled to bosonic thermal baths. This is so even 
when the bath interaction preserves the symmetry that 
protects the existence of edge states. As a consequence, 
these edge states decay in time into bulk states of a nor- 
mal insulator. Thus, a very fundamental question arises: 
is it possible to have stable topological insulating states in 
the presence of a thermal bath? The purpose of this let- 
ter is to explore this possibility by extending the concept 
of TI [2]-[5] to dissipative systems. Since for dissipative 
systems quantum states are generally mixed and charac- 
terized by a density matrix operator p, we shall refer to 
these as density matrix TI. 

For usual TIs, the TKNN invariant [6 provides a char- 
acterization of fermionic topological order. This is done 
in such a way that the transverse conductivity is writ- 
ten in terms of a topological invariant, the Chern num- 
ber, which may be related to an adiabatic change of the 
Hamiltonian in momentum space [7 . However the exten- 
sion of this invariant to density matrices is not straight- 
forward. Actually, the problem of generalizing geometric 
concepts as distances or geometric phases to generally 
mixed states is highly non-trivial ^8tiT2] . We address this 
problem and construct an observable that detects the 
topological order of a TI even if it is not in a pure, but 
in a general quantum mixed state. Moreover, when this 
general quantum state is of the form of a Gibbs state, 
we study the relation between this topological observ- 
able and the conductivity, and show that it reduces to 
the usual notion of TI in the limit of low temperature. 
However, we stress that the notion of a density matrix 
TI is far more general, as we shall see. 

This letter is completed with a supplementary infor- 
mation document (SM) where we include most of the 
technical details. 

Chern Connections for Density Matrices. — The physi- 
cal problem is defined as follows. Let be the sys- 




y 



Figure 1: Pictorial image of the action of a band Liouvillian 
C = Ck- The vertical lines denote the only possible pro- 
cesses involving the (initially empty) conduction band and 
(initially filled) valence bands, i.e. the ones where the mo- 
mentum k is preserved. The violet fog represents some bath 
at a certain temperature T which mediate such a processes 
(see Eq. ^) and the plane indicate the Fermi energy Ep- 



tem Hamiltonian representing a certain TI. This could 
be constructed in an arbitrary spatial dimension, but we 
shall restrict in what follows to the class of time rever- 
sal broken (TRB) insulators in two spatial dimensions. 
Furthermore, the TI will be subjected to the action of 
dissipative effects due to a thermal bath represented by 
a Hamiltonian i^b- This bath could be general enough so 
as to comprise fermionic or bosonic degrees of freedom 
and we assume it is initially in a thermal or Gibbs state 
at a certain temperature T. The system-bath interaction 
is described by the Hamiltonian i^s-b- 

We consider that the state ps of the TI undergoes 
a time evolution satisfying some Lindblad dynamical 
equation [13nl6j (unless otherwise stated, natural units 



h = ks = I are taken throughout the paper): 



dt 



(1) 



where C is the so-cahed Liouvillian operator, which is 
composed by a first term representing the Hamiltonian 
evolution in the absence of system-bath interaction and a 
second term accounting for the effect of the bath dissipa- 
tion. Concretely, we shall assume that C is of the Davies 
type, obtained under the assumption of weak system- 
bath coupling [17 . 

We are interested in searching for sufficient conditions 
that the Liouvillian dynamics ([T]) must satisfy in order to 
preserve the TI phase. In the absence of dissipation, we 
know that a key ingredient is that the TI Hamiltonian 
is a band Hamiltonian that satisfies the Bloch theorem 
and can be decomposed as = X^^gb z ^s(^) where 
k denotes a crystalline momentum. Thus, it is natural 
to restrict our attention to Liouvillian evolutions satisfy- 
ing a similar condition, jC = J^fceB z ^^^5 shall refer 
to these as band Liouvillians. Basically, a Liouvillian of 
this kind describes processes in such a way that the mo- 
menta of the fermions are not changed (up to a vector 
G of the reciprocal lattice), a pictorial image is sketched 
in Fig. [l] As a consequence, they present invariance un- 
der space translations, which imposes the following con- 
straints (see section [T] of SM) on its associated steady 
states pssXipss) = : 



Ao|0)(0|- 



a, 13 



la,fc)(l/3,fc|- (2) 



Here, \la,k) = \^cx,k) denotes a particle of type a in the 
Bloch state with momentum k associated with the sys- 
tem Hamiltonian i^g, where account for other in- 
ternal indexes, such as band indexes (like valence and 
conduction band) or spin indexes. The coefficients A^^ 
depend on the particular steady state as a result of the 
dissipative dynamics ([T]). For instance, if the steady state 
turns out to be a thermal state density matrix, then they 
will be given by Gibbs weights ([8|. 

In order to construct a topological invariant for the 
generally mixed state p^g, we cannot use the usual Berry 
connection as in the formulation of the Chern number, 
because it is defined just for pure states. Regarding den- 
sity matrices, there is not a unique natural extension of 
the Berry connection and the Berry phase [8--,-12^. How- 
ever, the band Liouvillian structure allows us for the con- 
struction of a Berry-type connection for the density 
matrix steady states, in such a way that the integral of its 
curvature form, F^-^ gives a topological invariant which 
we refer to as density matrix Chern value. Thus, one of 
the main results of this work is the construction of this 
invariant n^j^, which characterizes the topological struc- 
ture of insulators in the presence of dissipation and can 




Figure 2: Colormap depicting the Chern value in the Haldane 
model with dissipation, for different values of 0, M and bath 
temperature T (in units of t2 — 1). As T increases the Chern 
value decreases (in absolute value), and for T = we recover 
the phase diagram obtained by Haldane [18]. The dashed 
black lines enclosures the region displaying topological order 
at T = 0, so that all nonvanishing Chern values are inside 
of this region for any T. Approximately T — 1 and T = 5 
correspond to less than 10% and 50% of the gap respectively. 



be written as 

Oi 

+h^L ^^^'^"^ ^^^^^ " ^^'^"^ ^"^^^^ 

(3) 

Here, are the eigenvalues of p^^ in ([2| and Af{k) are 
the Berry connections of each band. We refer to the 
SM (section [2| for the detailed derivation and notation. 
A non- varnishing witnesses a topological non-trivial 
order present in p^g. When the steady state is a pure 
Bloch state p^g = |'WQ;,k)('^a,k|, we recover the standard 
TKNN topological invariant (Chern number). 

The density matrix Chern value n^^^, Eq. ([3|, has two 
different terms. The first one is a weighted integration 
of curvatures for different bands. This term has no topo- 
logical meaning on its own and it does not distinguish 
between phases with or without topological order. The 
second term represents a correction to the value given by 
the first one that provides the topological character to 
Hqyi- addition, both terms have a physical meaning 
which will be explained later on. 

The name Chern value responds the fact that despite 
its topological origin, it may not be an integer for a 
general mixed state. The reason is very fundamental; 
the space of density matrices p is a convex space, which 
means that a convex combination of density matrices pi 
and p2, p = pipi +P2P2 is also a mixed state. Due to 
the abelian character of the curvature form Ff- (see sec- 
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tion[2]of SM), n^j^ = pin^\^ P2^c\.' Therefore, since 
the weights pi,P2 ^ with pi + P2 = 1, then n^j^ G IR 
as weh. Nevertheless, this will be not an obstacle to use 
the Chern value to detect topological properties in in- 
sulator states. We will see this explicitly in the case of 
the Haldane model in 2D which is a prototype of TRB 
topological insulator ^18j . 

Band Liouvillian for the Haldane Model. — We can apply 
our previous formalism to the Haldane model of 2D TI. 
This is a graphene-like model based on a honeycomb lat- 
tice with nearest-neighbors and next-nearest-neighbors 
couplings. For periodic boundary conditions the Haldane 
Hamiltonian in the reciprocal space is given by 

Hs= Yl = E E^cick+Eldidk. 

keB.Z. ^ ^ keB.Z. 

(4) 

Here, and bk correspond to the two species of fermions 
associated with the triangular sublattices of a honeycomb 
lattice, and Ck and dk are the femionic modes which di- 
agonalize the Hamiltonian with eigenvalues and E2 
respectively. For more details we refer to section [3] of the 
SM. 

We shall assume a local fermionic bath model, with 
quadratic coupling of the form 

i^s-b := ^^g'ial^Al^ar^Al^^bl^Bl^br^Bi!), (5) 



where r denotes the point in the sublattices and 
and B^ the bath fermion operators coupled with the two 
species ar and br respectively. 

The detailed derivation of the Liouvillian equation 
(master equation) in the weak coupling limit for this sys- 
tems is explained in section[4]of SM; the final result turns 
out to be 



dt 



^Ck[Ps{t)]=^(-i[Hk,Ps{t)] 



+ 7(i?f)nF(i?f)(4ps(t)cfc - lickclpsit)} 

+ ^{E^)np{E^)(dip,{t)dk - \{dkdlps{t)}) + 
+ ^{E^)[1 - nFiE^)](dkPs{t)di - l{didk,Psm) 



(6) 



Here, 



e/^^ + 1 



, j{uj) :=27r J{uj), (7) 



where J{uj) is the bath spectral density. 

It is important to emphasize that this Liouvillian 
(|6| fulfills the conditions of a band Liouvillian C = 



J2keBZ ^f^' Moreover, it is quadratic in fermionic op- 
erators and its unique steady state is the Gibbs state 
{f3 = 1/T) 



-/3H, 



1 + e-^^i 



-(3E^dldk 



that has the form of pi corresponding to a band Liou- 
villian. Note that in the limit T ^ 0, p/3 approaches the 
Fermi sea where the lower band (c^) is fully occupied and 
the upper band {dk) is completely empty (see section [s] 
of SM for more details). 

Chern Value of the Steady State. — We have obtained 
that the steady state of the Liuvillian ([6| is a product 
of states with well-defined momentum. Thus, the 
(parallel) transport along k -space of each of these states 
is well defined, and hence, the state characterization by 
a density matrix Chern value (|3| is possible. 

For the sake of computation, note that p^^ 
is diagonal in the occupation basis p^^ = 
En,mem}Pnm\^^^)k{m,n\, whcrc 

|00)fc = |0)|0), \io)k = \uc,km, 

|01)fc = \0)\ud^k), = \uc,k)\ud,k)' 

The vacuum |0) has no particles and does not depend on 
k. If we define the geometric connections for the lower 
(c) and upper (d) bands: 



Af{k) := i{ua,k\diUa,k), a = c,d, 



(9) 



then it is possible to express the connection A^{k) in 
terms of the previous ones (see section [s] of SM): 

(it) = nF{E^)A'^{k) + nF{E^)Af{k). (10) 

Note that in the T ^ limit, we recover the standard 
Berry connection A^{k) (^), as the steady state ([8| 

approaches the Fermi sea with the fully occupied lower 
band. 

The Chern value can be now computed by integrating 
the curvature form of A^{k) (or by using the simplified 
expression ([3|). The colormap in Fig. [2] represents the 
Chern value for different values of M and (j) and different 
bath temperatures. Note its nice properties: it is zero 
for any choice of M and (j) for all temperatures if it is 
zero at T = 0. This manifests that the topological order 
cannot be created by increasing temperature. Moreover, 
as T increases, the absolute value of the Chern value 
decreases, and in the limit of T ^ 00 we obtain 
for all M and (j). This is in agreement with the common 
intuition that at infinite temperature any kind of order 
should be spoiled. 

Dissipative Edge States and Master Equation. — A phys- 
ical signature of a TI phase is the existence of gapless 
(metallic) edge states. Thus, once we have mathemat- 
ically characterized the phase diagram of the Haldane 
model under dissipation by means of the density matrix 
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Figure 3: Occupation along the direction ai for all particles 
with momentum k2 — fy^? M = and — 7r/2), for dif- 
ferent values of T (in units of ^2 = !)• Note the presence of 
edge states at finite temperature, Eq. ( p^ , in the positions 
1 and 50 along the direction ai. Note the presence of edge 
states at finite temperature, . However, as the temperature 
T increases significantly the population of the edge modes 
becomes similar to the population of the bulk modes. 



Chern value, we wonder about the fate of the chiral edge 
states of the Haldane model at finite temperature. 

To that aim, we consider the Haldane model placed on 
a cylindrical geometry, where we take periodic boundary 
conditions just along one spatial dimension, say a2. In 
such a case the momentum k2 along the a2 direction is 
a good quantum number and the Haldane Hamiltonian 
can be diagonalized obtaining (see section [s] of SM for 
more details) 

H,= Yl ^(^2)= E E'-fl,J^,k,, (11) 

/C2GB.Z. /C2GB.Z. 

Here, the diagonal modes fmM ™x both species of 
fermions am,k2 and bmM- 

By imposing this geometry also in the interaction 
Hamiltonian ([5| , we derive the following dynamical equa- 
tion for the system (see section [4] of SM) 

^= cMt)]= E (-mk2),Ps{t)] 

' \ '> {m,k2) 

m 



where 



VK[ps{t)] := Kps{t)K^ - ^{K^K,ps{t)}. 



(12) 



(13) 



Again, the Gibbs state at the same temperature as the 
bath is the unique steady state of Eq. (12), 



P/3 



-/3//(/C2 
Zk2 



(14) 



with Zk^ = Tr [e-/^^(^2)] ^ Therefore, as long as the val- 
ues of M, t2 and (j) are such that the system exhibits 
topological order (see Fig. [2|, two of the modes which di- 
agonalize each H{k2)^ say /(l,/c2) and f(R^k2)^ correspond 
to edge states and the Gibbs state is a tensor product in 
/c2 of states of the form 



where 



-pH{k2) 



P^(^2)0P^^^^^2)0P^(^2), (15) 



^ ^ Q-I^^L,R{k2)fl^^^^^^^f{^L,R,k2) 

are Gibbs states of the edge modes. However as temper- 
ature increases (see again Fig. [2| t he edge modes in the 
steady state of the Liouvillian ^12| become delocalized. 
Then the components p^'^(A:2) approach to a maximally 
mixed state with vacuum |0). In such a situation, there 
is not a way to distinguish the edge from the bulk modes, 
because for all k2 the occupation along the direction a\ 
becomes the same and equal to 1/2. We illustrate this 
behavior in Fig. [3] 

It is worth stressing that in order to study the dissi- 
pative effects on the cylindrical Haldane system we must 
determine how the boundary conditions of the system 
affect the dissipator operator. To that aim we need to 
specify how the dissipator is generated (constructed). In 
our case this is the result of the weak interaction of the 
system with local baths. However there are also pos- 
sible scenarios where the dissipator is the effective re- 
sult of other external interactions (see for example [19]- 
23 ). Since the process of "opening" or "closing" a system 
has a physical meaning, we stress that we need to know 
how the dissipator is physically generated to obtain its 
"open" and/or "closed" counterpart. Note that a dis- 
sipator with periodic boundary conditions, may split in 
different and nonequivalent dissipators once the system is 
opened along some direction if it is generated in different 
ways. 

Quantum Hall Conductivity and Chern Value at Finite 
Temperature. — We can obtain further physical meaning 
and implications for the density matrix Chern value (|3| 
by studying the (quantum Hall) transverse conductivity 
axy and its relation to the thermal edge states obtained 
for the Haldane model. Using the Kubo formula [24^ in 
linear response theory, it is possible to derive an expres- 
sion for the transverse Hall conductivity at finite temper- 
ature [25l: 



27Th ^ 



f nF{El)F^y{k)d^k. (17) 
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the conductivity: 



Figure 4: Colormap depicting the conductivity Eq. (17) for 



different values of 0, M and the bath temperature T (in units 
of t2 = 1). As T increases the conductivity decreases (in 
absolute value), and for T = we recover the Chern number 
result, Fig. |2] The dashed black lines enclosures the region 
displaying topological order at T = 0. Thus, contrarily to the 
Chern value, the conductivity is not topologically invariant 
for finite T, as it does not vanish for every point outside this 
region for any T. 



Note that this expression [26lis different from the one 
obtained for the Chern value pi). Indeed, the conduc- 
tivity is not topological at finite temperature, as shown 
in Fig. |4] where non-zero Hall conductivity appears in 
regions outside the topological regime, in contrast with 
Fig. |2j Nonetheless, both quantities can be related by 
means of the following equation: 



— n^ 



27rh 



(18) 



d'k. 



The second term on the right hand side of (18) is the 



same one that appears for the transverse conductivity of 
a normal insulator with an applied magnetic field (or 
a pseudo-magnetic field as for the Haldane model) at 
T 7^ 0. It corresponds to the conduction by thermal 
activation of excited electrons in the bulk. For instance, 
for parameters ti = 4, ^2 = 1,^ = f and M = 6 in 
the Haldane model, this term is the only non-zero contri- 
bution to the conductivity, as the system is outside the 
topological regime, and so = 0. 

Notwithstanding, the first term on the right hand side 
of (18), which is nothing but the Chern value previously 



defined, represents a contribution due to the topological 
nature of our system and the presence of conducting edge 
states. For parameters ti = 4, ^2 = 1, ^ = f and M = 
in the Haldane model, the system is within the topolog- 
ical regime and this new term shows up. Note that at 
T = we recover the well known TKNN expression for 



>o h 



(19) 



where uch denotes the standard Chern number. 

Conclusions. — We have studied topological insulating 
phases in the presence of dissipation. After introducing 
the notion of band Liouvillian, we address the charac- 
terization of the topological order of its steady states by 
resorting to the density matrix Chern value, a topolog- 
ical invariant that is an extension of the Chern number 
for pure states. The Haldane model of a 2D TI in con- 
tact with a thermal bath offers a nice testbed to study 
these phenomena. More concretely, we compute phase di- 
agrams at finite temperature based on the Chern value, 
and corroborate that topological order decreases as the 
bath temperature increases. Thus, from a topologically 
disordered state it is not possible induce a topologically 
ordered phase just by warming the system. However a 
topologically ordered state may remain ordered at finite 
temperature T except at the limit T ^ oo. This has 
to be compared with the previous study pi where the 
topological order turned out to be lost for a dissipative 
system in the presence of noise not generated by a band 
Liouvillian. Our results may also have direct applica- 
tion in recent studies regarding dissipation on Majorana 
fermions in topological superconductors [20^, ^27^ . 

Complementarily, we study the properties of the Hal- 
dane model coupled to a thermal bath under cylindrical 
boundary conditions. We find that the the steady state 
splits in three different, generally mixed substates Eq. 
(15), two of them are associated with creation or annihi- 



lation of fermionic gapless edge modes, and the other one 
accounts for the same process just in the bulk modes. In 
the limit T ^ we recover the properties of the usual 
Haldane model. 

Finally, we examine the relation between the density 
matrix Chern value and the conductivity at finite tem- 
perature. We show that the latter is not topologically 
invariant, in contrast to the Chern value. This fact is 
due to the presence of an extra term which accounts for 
the conductivity generated by thermal activation of elec- 
trons in the bulk, which has not a particular topologi- 
cal meaning and is present in normal insulators. In this 
regard, provided that the gap between conduction and 
valence bands is large enough, the Chern value may be 
approximately estimated by measuring the conductivity, 
as the thermal activation would be very small. However, 
recent results [28 , suggest that a direct measurement of 
the density matrix Chern value could be possible in op- 
tical lattice realizations. 
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SUPLEMENTARY MATERIAL 

1. Band Liouvillian Dynamics 

There are some natural requirements to be imposed on 
a Liouvillian C in order to generalize the theory of TIs as 
developed for band Hamiltonians = X^fceB z ^s(^) 
where k denotes a crystalline momentum. In analogy 
to these, we introduce the class of band Liouvillians, 
which by definition are those that can be decomposed as 
jC = X^fcGB z ^^8' where each Cj^ only involves fermionic 
operators with the crystalline momentum k. A band Li- 
ouvillian presents translational invariance, as every Ck is 
translationally invariant: 

T{a)Ck{p)T\a) = Ck [T{a)pT\a)] , (SI) 



where T{r) 



-iak 



is the operator that translates a 



point with coordinate r to the point r + a on the lattice. 

An analogous to the Bloch theorem for this kind of 
Liovillians characterizes its steady states. 

Theorem 1 Consider a band Liouvillian C = 
^keB z ^k- U ecich Ck has a unique stationary 
state, it has the form 



Ao 



E 

k,a,(3 



A^;3|la,fc)(l/3,fc| 



(S2) 



Here, a and j3 denote additional quantum numbers (band 
indexes, spin indexes, lattice indexes, etc.), and \la,k) = 
\ua,k) denotes a particle in the Bloch state with momen- 
tum k and additional quantum number a, T{r)\ua,k) = 
e'^'^'^lua^k)' 

Proof: Consider p^s to be the steady state of some Ck- 



Ck{Pss) = 0. 



(S3) 



By applying the translation operator on both sides and 
using (SI), we obtain 



Ck [T{a)pssTHa)] = 0. 



(S4) 



Thus, := T{a)pssT\a) is also a steady state of the 
system. Since by assumption pss is unique, = pss, 
so that [T(a), pss] = and T{a) and p^s share the same 
eigenvectors. □ 
It is worth noticing that as a difference with the case 
of pure states, the translational symmetry of a Liovil- 
lian does not necessarily implies steady states with well- 
defined crystalline momentum k. They can be a con- 
vex mixture of states with different well-defined momenta 
(S2). However, since the subspace with well-defined mo- 



mentum k is invariant under the action of Ck^ if the ini- 
tial state has a well-defined momentum k (for instance a 
particle with well-defined momentum in one of the bands 
of the Hamiltonian), then the steady state under Ck will 
have well-defined momentum as well. 



Ao|0)(0| 
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(S5) 



Thus, the steady state of the total Liouvillian C = 
XlfcGB z form of 



(S6) 



Purified Topological Invariants for Band 
Liouvillians 



Purification is a method that allows us to extend quan- 
tities defined for pure states to general mixed states. For 
a density matrix p acting in a Hilbert space 'H, a purifi- 
cation |<I>^) is a pure state in an extended Hilbert space 
eUA^U such that 



p = TTA{m{^' 



(S7) 



In other words, mixed states can always be seen as pure 
states of a larger system such that we only have access 
to partial information of it. 

Given some p there are infinitely many states |<I>^) 
which fulfill (S7). Without loss of (mathematical) gen- 
erality, it can be assumed that the ancillary space 1-La 
(which physically may represent the thermal bath) has 
the same dimension d as 1-L for the system, then any pu- 
rification 1$^) can be written as 



\<^>p) = {UA^p)\n), 

where Ua is a unitary operator, pp^ = p and 

d 



(S8) 



(S9) 



is a maximally entangled state, with an orthonor- 

mal basis of H. From the Schmidt decomposition of |<I>^) 
it follows that ( S8 ) is the most general form for a purifi- 
cation of p [1 . 

By using the spectral decomposition of p = 
we may write p as 



(SIO) 



where {Iv^a)} is also an orthonormal basis which is con- 
sidered to be arbitrary. Therefore, given some p there is 
a freedom for the choice of Ua and the basis {l^/^^a)} for 
its purification \^^). 

The next theorem is the main result of this section. 

Theorem 2 Consider the steady state of a band Liouvil- 
lian, Eq. (S6); we define a Berry-like connection for p^g 



through one of its purifications |$^) as 

Al{k) := liKmi), for pi = TVA ■ 

(sii) 

Here, the notation is di := dk^- Under the assumption 
that Ua ctnd {\^i)} are independent of momentum k of 
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the steady state, the connection (Sll) is unique and does 



not depend on the purification. Explicitly, it takes the 
following form in terms of the spectral decomposition of 

Pss = 12o.Pa\i'a,k){i>a,k\: 

A^{k) = iY,p1{^o.,k\di^o.,k)- (S12) 



Proof: Indeed, the general form ( S8 ) for a purification of 
reads 



where we h ave used (SIO). Taking the derivative di := 
in (|S13|) and computing the overlap 



a, (3 



9i\ Pa ] (1 ^ \^l3){i^^,k\^a,k){(Pa\) 



\n). 

(S14) 



Since ^ A)\n) = Tr(A), we obtain 

^fW =iE^(^^^) ^Pa{^a,kma,k). (S15) 

Q! ^ ^ 

which is independent of Ua and Moreover, note 

that in the pure state case p^g = we recover the 

Berry connection A^{k) = Ai{k) = i{^pk\di^lJk) [2^- In 

addition, the first term on the right hand side of (S15) 
vanishes by taking into acco unt t hat ^^p^ = 1- Hence, 
(S15) can written simply as (S12). □ 



Once this purified connection (S12) is defined, we may 



obtain the (abelian) curvature form through 



(S16) 



and construct a density matrix topological invariant 
associated with the steady state p^g via the first Chern 
class |3] of this connection. 



271 ' 



-Tr 



T2 



F[.{k)dki Adkj 



(S17) 



It is convenient to compute the different contributions 
that appear in the exphcit expression of (S16) after using 
(ISni). We obtain 



E 



[p'^Ftj + {d,p'^) A^{k) - AUk)] . 

(S18) 

Note that from this equation is not manifestly clear the 
U{1) gauge invariance of the curvature, but this can be 
proven by performing a gauge transformation and mak- 
ing use of the property X^q^^S ~ ^- However, the gauge 
invariance is more straightforwardly shown by n oticing 
the clear gauge character of the connection (Sll). 




Figure 81: System of coordinates {a i, 02} taken to write the 
Haldane Hamiltonian in real space (S19) . The solid white 



and black circles denote fermions a and h respectively, and 
the green enclosures highlight the two-site unit cell. 



3. Haldane Model Geometries 

In order to write explicitly the Haldane Hamiltonian 
j4j in real space we will consider the system of coordinates 
{ai,a2} represented in Fig. SI, with ai = ^(3, >/3) and 



a2 = ^(— 3,>/3) for lattice spacing a = 1. We write 
<^(m,n) (^(m,n)) fo^" ^hc fcrmiouic operator of kind a (b) 
in the position = mai + na2. Then, the Haldane 

Hamiltonian in real space reads 



(^'2 \-'^lm,n)'^{m,n) ^(m,n)^(™,n)] 
+ *1 ['4m,n)^(rn,n) + «(m+i,„)^m,n) + a(™,„)^(m,n+l) 



t2 



e"^«(m,n)«("»+l.n+l) + ^ "^«(m,„)«(m+l,n) 



+ e "^a(™,„)«(m,n+l) +e "^&[„_„)&(m+l,n+l) 



h.c. . 



(S19) 



a. Toroidal Geometry 

By taking periodic boundary conditions in both spatial 
directions, we may write the Hamiltonian (S19) in the 



reciprocal space using the Fourier-transformed operators 
1 J2 e^^"-(-")afc, (S20) 



keB.Z. 



fcGB.Z. 



(S21) 
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where B.Z. stands for Brillouin zone which is a hexagon 
with vertices in the k = (/ci, /C2) points 

(0'^)'(¥>^)-(f--^)- (S22) 

(0.-3^)'(-¥.-^).(-¥'^). 

and N is the total number of two-site unit cehs. Thus 
the Haldane Hamiltonian is rewritten as 



(S24) 



Here, 



Hii{k) = M + 2t2^cos[0+(/c-6i)], 

i 

i 

H22{k) = -M + 2t2^cos[</.- (fe-bi)] 



(S25) 



with 



61 = -(3,v^)/2, 

62 = (3,-\/3)/2, 

63 = (0,^/3). 

By diagonalizing the matrix H{k) we obtain 

where the eigenvalues are given by 

E^^ = 2t2(cos(/)) Si{k) T ^/^k) 

with 



(S26) 



(S27) 



(S28) 



A{k) := ^tl[3^2di{k)] 

+ 4t^(sinV) Mk)]^ - ^Mt2{sm(^)62{k), (S29) 



and 



Si (k) := cos k -bi = 2 cos ^ cos 



\/3fe2 



cos 



^2(^) := EiSin/c-6i = -2 cos ^ sin ^ + sin V^A:2. 

(S30) 

The operators and are related to and via the 
canonical transformation 



Vi + l^fcP v-1 ^iJ 



where 



Xk := 



kai 



M - 2t2(sin0) S2{k) + v'A(fe) 



(S31) 



(S32) 




Figure S2: Energy bands of the Hamiltonian (S35) as a func- 
tion of the momentum /c2. The red lines correspond to the 
edge state modes. The rest of parameters are M = 0, (/) = 7r/2 
and ti — A (in units of 12 — 1). 



h. Cylindrical Geometry 

In this case we take periodic boundary conditions along 
the direction a2 and open boundaries along ai, so that 
we work with a cylindrical configuration. The inverse 
Fourier transform along the a2 direction of the fermionic 
operators is given by 



^(m,n) 
^(m,n) 



k2n\a2\ 



2 /C2GB.Z. 



(^{mM)-> (S33) 



2 /esGB.Z. 



where the Brillouin zone corresponds to the interval 
k2 e (-7r/|a2|,7r/|a2|) = (-tt/VS, tt/V^), and N2 is the 
number of two-site basic cells along the direction a2. By 
using these equations in the Hamiltonian (|S19[) we obtain 



Hs= E([f+*2 cos (73^2 



^(m,/C2)^(^'^2) 



[-f +t2 cos(V^A:2 + 0) 



(^/3fe2+(/>) 



^(m,/C2)^(^+l'^2) 



y ^(m,fe2)^(^+l'^2 



h.C. 



(S35) 



This Hamiltonian has the structure = 
X1/C2GB z ^(^2), therefore we can diagonalize by 



diagonalizing each H{k2). In Fig. S2 we have depicted 



10 



the behaviour of the eigenvalues of H{k2) as a function 
of k2- The red hues connecting the upper and lower 
bands correspond to the edge state modes, which are 
localized on the edges of the direction ai. 



4. Derivation of the Master Equation for the 
Haldane Model 



In this section we derive the dynamical equation (mas- 
ter equation) for the Haldane model coupled to a thermal 
bath. We assume the usual condition of weak system- 
bath coupling, which is a standard assumption for ther- 
malization. 

The total Hamiltonian of the problem considered reads 
as follows 



where Trb denotes the trace over the bath degrees of 
freedom and is the initial state of the bath, which we 
assumed to be the Gibbs state 



(S42) 



a. Toroidal geometry 

We consider periodic boundary conditions and take 
Fourier transforms in the Hamiltonian (S36). For the 



first term we obtain ( |S24[ ), for the second we have 

Hy, = Y,e'(^^A'^ + B^Bi), (S43) 



H:=Hs + Hk + Hs-h- (S36) 

The first term, i^g, is th e Haldane Hamiltonian ( |S19[ ). 
The second term in (S36), i^b, is the free Hamiltonian of 
the local baths. 



(S37) 



where A and B stand for independent fermionic bath op- 
erators that satisfy the canonical anti-commutation rela- 
tions {Al,A^\} = 6n,n'^ij , {Al,A{,} = and analo- 
gously for Bl^. The index r denotes the position of the 
local bath on the lattice and i runs over the bath degrees 
of freedom. Also, represents the energy of each mode 
i of the bath which is assumed to be independent of the 
lattice site. Finally, the third term in (S36), i^s-b, is 
given by 



Zs-b • = 



g\al (g) a; + (g) AjJ -^bl^B^^br^ B'J), 

(S38) 

and describes an exchange of fermions between system 
and bath mediated by a coupling constant which may 
depend on the specific mode i of the baths. 

The total dynamics of system and bath is given by the 
Liouville-Von-Neumann equation: 



dp 
'dt 



i[H,p], 



(S39) 



After taking the interaction picture with respect to Hq = 

^ = -i[i7,_b,p] with |^^_^^^,i.o.^^_^e--o.. 

(S40) 

For small i^s-b the system dynamics is approximately 
given (see [5]-[8]) by the equation 



dps 
dt 



■ / ds Trb[ils-b(i), [Hs-h{t - s),~Ps{t)®pi]], 
Jo 

(S41) 



and finally for the interaction term 



s-b 



i,k 



g\al (g) + afc (g) ^hl^Bl^hk^Bl 



(S44) 

Let us stress again that the strength of the coupling to 
each mode of the bath is represented by g^ . This, anal- 
ogously to the energy of each mode e*, is taken to be 
independent of the lattice site and of the type of bath A 
or 5, which is rather natural. 

In terms of the operators Ch and dh the Hamiltonian 
(IS44|) reads 



i^s-b = ^9\cl^Ci^Ck®Ci^ ^dl^Di^dk^D^). 

(S45) 



i,k 



where 



\Xk\ 



1 Xk 



(S46) 



are new fermionic modes of the bath. Moreover note that 



i,k 



(S47) 



Now, it is easy to write i?s-b in the interaction picture 
and apply the formula (S41), which can be quite simpli- 
fied by using that 

TMCilcipp) = Tr,,{Di\Dipp) = nF{e%j5k,k', 
TrkiC'llDlpp) = Tv^iDilClpp) = 0. 



Here, np^E) := 



Q^Ej^-^ is the mean number of particles 
of the Fermi-Dirac distribution, where we have taken the 
chemical potential p to be at the origin of the energy. In 
the continuous limit for the baths degrees of freedom we 
have 



^ j deJ{e)f{e), 



(S48) 
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for any function /(e), where J {00) is the so-cahed spectral 
density of the bath. Thus, the Sokhotsky's identity 



^ dre^^^=7rJ(cj)+iPv(^^ 



(S49) 



ahows us to simphfy further the final expression, which 
after a bit long but straightforward computation reads 



dt 

+ 7(^f )[1 - nF{E^)] (c^ps W4 - Psd)}) 
^^{E^)nF{E^)(dip,{t)du - \{dudlp,{t)}Y 

+ 7(^2)[1 - nF{E^)](duPs{t)di - \{didu. p,{t)} 

(S50) 

in Schodinger picture, where 7(0;) := 2ti J {uo). In addi- 
tion, in this equation we have neglected the imaginary 
parts of Eq. ( |S49[ ) because they represent just a small 
shift of energies which does not affect to the dissipative 
process [9]. 



h. Cylindrical geometry 

In this case, we take Fourier transform along the direc- 
tion a2 in (S36). Thus, the Haldane Hamiltonian reads 



as ( |S35[ ), the bath Hamiltonian becomes 

(S51) 



/C2 i,rn 

and the interaction Hamiltonian 



The diagonal modes f{mM) ^(^2) = 

T.m^rn{k2)fl^^f^^^f(mM)^ whcrc £^^(^2) is depicted 



in Fig. S2, are related to C(^rnM) some unitary 
transformation 

C{mM) =^^m,if{iM)^ with ^W^^^W^^^, = Sm,m' - 

^ (S54) 
By using this equation in ( S53 ) and after a bit of algebra 
we arrive at 



(S55) 



i,m,/e2 



where 



(S56) 



are new fermionic bath modes. 

Following the same steps as for the toroidal geometry, 
it is not difficult to obtain the master equation of the 
cylindrical array. 



dt 



fe2GB.Z. fe2GB.Z. ^ 

' \ '> {m,k2) 

m 



(S57) 



where 



VK[ps{t)] ■■= Kps{t)K^ - ^{K^K,ps{t)}. (S58) 



5. Steady State for the Haldane Model 



a. Toroidal Geometry 



/C2GB.Z. i,m 



The steady state of the previous band Liouvillian ( S50 ) 
is the Gibbs state of the Hamiltonian i^g, 



(S52) 



where m runs from 1 to the number of two-sites basic 
cells along the direction ai, Ni. 

We may collect the operators and bm of each site 
in a new operator Cm where Ci := ai, C2 := C3 := a2, 
C4 := 62 and so on. The same can be done for the baths 
operators A and B with the notation Cm- Then, the 
interaction Hamiltonian ( S52 ) is written as 



(S59) 



/C2GB.Z. i,m 

where now m runs from 1 to 2Ni. 



^«-b= Yl IZ^'Hm,fe2)^^(rn,fc2)+^(-,fc2)^^(l,fc,) 

(S53) 



where (3 = 1/T, with T the temperature of the fermionic 
bath, and Z = Tr(e~^^^) the partition function. To 
prove this, first note that [Hs^ p^] = so we just need 
to care about the dissipator in ( |S5Q| ). In addition, since 
the number operators c^c^ and d]^dk commute with c^,, 
Cfc/, (i^,, and d^^ if k ^ , and we are left only with the 
part where the crystalline momenta in and p^^ are the 
same, as the others trivially vanish. Since 



e^^nF{E) = [l-nF{E)], 



(S60) 
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and 



it is easy to proof that 



-(3E^dldu 



(S61) 



CkiPss) = fiF{E)(ckpici - l{cick,pi} 



[1 - nriE)] (ckpici - ^{^c^, p^} 



nF{E)[dkpidi - l{didk,pi} 



[1 - n^m (^dj^pidi - ^ididj^, pi}^ = 0. 

(S62) 



Moreover, the state p^ is the uniq ue st eady state of ( |S5Q| ) 
as the interaction Hamiltonian (S38) satisfies the irre- 



ducibihty condition presented in jlO . 

In order to analyze some properties of let us write 
the density matrix p^^ in the occupation basis of the two 
bands for a fixed momentum k. The basis reads as \ij)k^ 
where z = 0, 1 and j = 0, 1 stand for the occupation of 
one-particle state in the lower and upper band respec- 
tively. Thus, p^g is a 4 X 4 diagonal matrix. 



diag (poooo^PioiO'Poioi^Piiii) 



(S63) 



with 



Poooo 

rfolo 
Poioi 
Pun 



(1 



• Poooo^ 

• ^0000^ 



e-^^f)(l 



(S64) 



Since we set the origin of energy at ^0 = and also 
took the chemical potential p = in between the two 
bands, < 0, > 0- Thus, at the low temperature 
limit, the state p^g |10)fc as T ^ K. This means 
that at T = K the lower band is fully occupied and the 
upper band is completely empty, that is actually what 
one may expect. At the opposite limit, for T — > 00, the 
system gets completely mixed, as the four possible states 
can be equally populated by the environment. 



For the sake of clarity, we write the members of the 
occupation basis as 

|00)fc = |0)|0), 
|10)fc = \uc,k)\0), 
\Ol)k = \0)\ud,k). 

= \Uc,k)\Ud,k)' 



Then by using Eqs. ([9| and (S12), we obtain (note that 
the di\0) = as by definition the vacuum has no particles 
and so it does not depend on k) 

(S65) 

where the Berry connections for the lower c and upper d 
bands are provided by 

Af{k) = i{ua,k\diUa,k), a = c,d. (S66) 

Thus, the expressions for p^jj^i given in Eq. ( |S64[ ) lead to 

<(fe) = np{E'^)AUk) + np{E^)Af{k). (S67) 

b. Cylindrical Geometry 

Because of the same reasons as with the toroidal ge- 
ometry, the master equation on a cylindrical geometry 
(S57) has a unique steady state, which is the Gibbs state 



at the same temperature as the bath 



(S68) 



with Z/c2 = Tr [e /5^(^2)j^ Provided the system exhibits 
topological order, this state can be split as 



P^{k2) = ^-^ = p^(fc2)^P^"^^fc2)0p^(^2), (S69) 

where 



(S70) 



are Gibbs states involving just the gapless edge modes 
depicted in Fig. [82) 
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